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ABSTRACT 



Thanks to the ever increasing computational power and the development of more 
sophisticated algorithms, numerical TV-body simulations are now uncovering several 
phcnomenological relations between the physical properties of dark matter haloes in 
position and velocity space. It is the aim of the present work to investigate in detail 
the dynamical structure of dark matter haloes, as well as its possible dependence on 
mass and its evolution with redshift up to z = 5. We use high- resolution cosmological 
simulations of individual objects to compute the radially-averaged profiles of several 
quantities, scaled by the radius i? m ax at which the circular velocity attains its maxi- 
mum value, Vmax- We investigate the circular velocity profile, the dark matter density 
and its logarithmic slope, the position and velocity of the centre of mass on different 
scales, the radial infall around the halo, its spin parameter, the radial and tangential 
components of the velocity dispersion, and the coarse-grained phase-space density. It 
is found that all the physical properties considered display a similar structure when 
expressed in terms of i? max and V ma , x . No systematic dependence on mass or cosmic 
epoch are found within -R maX j and all the different radial profiles are well fit by simple 
analytical models. However, our results suggest that several properties are not 'uni- 
versal' outside this radius. In particular, dark matter haloes should not be assumed 
to be in equilibrium beyond -ff max , especially at high redshift, where significant infall 
velocities can be measured. We therefore conclude that the whole dynamical structure 
of haloes, rather than just their density profile, is 'universal' within i? max at least up 
to z = 5, and that some dynamical process yet to be identified is responsible for such 
universality, rather than the details of the merger history of the halo. Memory of the 
initial conditions is only retained in the outer regions, where there is significant scatter 
from object to object, as well as clear systematic trends with mass and time. Thus, 
we argue that i? max and V^ lax provide a well-motivated, model-independent choice to 
characterize dark matter haloes, especially in comparison with the so-called 'viriaP 
mass or radius. 
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1 INTRODUCTION 

The internal structure of dark matter haloes is one of the 
key predictions of any model of structure formation. Al- 
though the complexity of the non-linear phase of gravita- 
tional collapse has so far precluded an exact analytical treat- 
ment, much progress has been possible over the last three 
dec ades thanks to cosm ological numerical experiments (see 
e.g. iBertschingerll 19981 . for a historical overview) . 

Indeed, perhaps one of the most well-known results 
of A-body simulations is the so-called 'universality' of the 
dark matter density profile, p(r), of the simulated haloes 
jNavarro et a.1.1 fl997f ) . Although the asymptotic behaviour 
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towards the centre (see e.g. iMerritt et al.l 120061 . and refer- 
ences therein, for a recent discussion) or the outer parts 
jPrada et al.ll2006t ) are still open questions, there is general 
consensus in that the radial density profile of dark matter 
haloes can be well fit by a simple analytical function with 
very few free parameters. Many of the functional forms pro- 
posed in the literature can be cast in the form 

P{r) = x" [1+% /r )T (1) 

where the characteristic density and radius po and ro are 
usually regarded as free parameters, and different prescrip- 
tions exist for the exponents a, /3, and 7. It is important 
to note that po and ro have been found to be correlated, a 
phenomenon referred to as the 'mass-concentration relation' 
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it has been claimed that, for any given object, the product 
pprp stays approximately constant throughout its evolution 
l|Romano-Di'az et alj|2007f ). 

In velocity space, dark matter haloes also seem to 
display a remarka ble degree of universa l ity, w i th the angular 
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20011 : IRasia et all 12004 IDehnen fc McLaughlin! 120051 ) pro- 
files being well described by simple analytical fits at 

2 = 0. 

Of course, the density profile and the velocity struc- 
ture of the haloes are not independent, and different cor- 
relations have been observed. For instance, a linear re- 
lation between the anisotropy parameter (3 and the lo- 
cal slo pe of the density profile has recently been advo- 
cated (|Hansen fc Mooreil2006l: lHansen fc Staddl2006h . and 
it has been arg ued l|Tavlor fc Navarroll200lf ) that the coarse- 
grained phase-space density profile, q(r), of galaxy-sized 
haloes follows a power law 



q{r) = 



P{r) 
<T 3 (r) 



(2) 



wi th exponent a = 15 /8 = 1.875. This result was confirmed 
bv lRasia et al.l (|2004T ) on the scale of galaxy clusters, where 
the value a = 1.95 was found t o pro vide the best fit to 
numerical data. lAscasibar et al.l (|2004l ) showed that equa- 
tion ([2| actually holds at all scales, from galaxies to galaxy 
clusters, when expressed in units of 'virial' quantities, 



9(r) = K„ 3 



V: 



( Rvir ) 



(3) 



where p c denotes the critical density of the universe, M vlr 
and iivir are the virial mass and radius of the halo, V^ ir = 
GM vil /R viv , K = io 146±0 04 , and a = 1.90 ± 0.05. The 
scaling of q with halo mass or v elocity dispersion has also 
been studied by s everal authors ( Dalcanton fc Hoganll200~il ; 
iDave et al.ll200ll : IPeirani et alJl2006h . and the evolution of 
this quantity durin g the merging hist ory o f the halo has 
been addressed by IPeirani et all [120061, and Hoffma n etaD 
l|2007t ). Using the Jeans equation. pTavlor fc Navarro! ()200ll ) 
computed the density profile expected for systems in equi- 
librium that fulfill equation ([2]). Families of analytical so- 
lutions ha ve later been discussed in detail for both the 
isotro pic (|Hansen| 120041 : lAustin et al. 20051: iBarnes et al.l 
20061) and the aniso tropic case ( Dehnen fc McLaughlinll2005l : 



Bar nes et al.ll2007h . 



To summarize, there is plenty of evidence suggesting 
that the present-day dynamical structure of dark matter 
haloes may be quite insensitive to the details of their ac- 
cretion history. More precisely, the radially-averaged profiles 
of 

(i) density (or mass) 

(ii) angular momentum 

(iii) velocity dispersion (radial and tangential) 

(iv) coarse-grained phase-space density 



are similar for all haloes, and can be described in terms of 
a few free parameters. In this paper, we use a set of high- 
resolution cosmological N-body simulations to investigate 
the dependence of the dynamical structure on halo mass and 
epoch, up to 2 = 5. Apart from the quantities listed above, 
we also consider the radial infall velocity and the offset be- 
tween the centres of mass on different scales. It is shown that, 
at all redshifts, all the relevant properties of dark matter 
haloes can be expressed in terms of the characteristic radius 
i?max at which the circular velocity, V c 2 = GM/r, attains 
its maximum value, Vmax- These two parameters are model- 
independent, they can be easily measured for any dark mat- 
ter halo, and they provide a reasonable indication of the 
extent of the inner region where the halo is in equilibrium. 

The paper is structured as follows: our numerical ex- 
periments are briefly described in Section[2] Results are pre- 
sented in Section [3] and the main conclusions are summa- 
rized in Section [4] 



2 NUMERICAL EXPERIMENTS 

The present work is based on two sets of A-body simula- 
tions carried out wit h the Adaptive Refinement Tree code 
|Kravtsov et al.|[l997T ) . Both of them are random realizations 
of a ACDM cosmology with parameters Q m — 0.3, £7a = 0.7, 
h = 0.7, and as — 0.9, and the box size is L — 80 h -1 Mpc. 
In one case, we have selected eight cluster-sized haloes at 
z = to be resim ulated by the multiple-mass technique 
ijKlvpin et al ] |200ll ). while six galaxy-sized haloes have been 
selected in the other set. Particle masses are m p = 3.16 x 10 s 
and 5.0 x 10 6 h -1 Mg, respectively, and the force resolution 
reached (two times the size of the highest refinement level 
cell) was 2.4 and 0.6/i -1 kpc in each case. Mass and force 
resolution are an important concern in order to obtain reli- 
able measurements of the internal structure of dark matter 
haloes. At high redshift, mass resolution becomes especially 
critical, since even the most massive objects are composed 
by a small number of particles. 

The mass distribution of this sample of 14 simulate d 
haloes has been investigated in Ascasi bar et al.l l)2007h . 
where it was compared with the predictions of the spher- 
ical infall model. Here we extend the analysis to the internal 
structure of these objects in velocity space. It is interesting 
to stress that our haloes have been randomly selected, i.e. 
no effort has been made to select isolated nor dynamically 
relaxed systems, and therefore they represent an unbiased 
sample of the average cosmic population. 

For each halo, we compute the position of the centre of 
mass iteratively, starting with an initial guess and a sphere 
of 30 — 300 h~ x kpc radius. When convergence is reached, the 
radius of the sphere is reduced by ten per cent, and the pro- 
cess is repeated until it encloses less than 10 3 particles. This 
procedure is almost equivalent to selecting the minimum of 
the g ravitational pot ential as the origin of coordinates (see 
e.g. 



Neto et all 120071 ). Then, Knax and -Rm ax are obtained 



from the circular velocity profile around that point. Radial 
bins are chosen so that the enclosed mass increases by 5 per 
cent (or 100 particles, whichever is larger). The innermost 
bin in all our profiles encloses more than 100 particles. 

To locate the first peak in V c , we examine the profile 
in order of increasing radius until the circular velocity de- 
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Table 1. Radius -Rmax (in h 1 kpc) where the circular velocity V m ax (expressed in km s 1 ) reaches its maximum, and number of 
particles contained within -Rmax, for our sample of simulated dark matter haloes, at the three redshifts considered. 



creases below 97.7 per cent of the maximum value found so 
far; in other words, 



log(Vc) =0.99 log( Vmax). 



(4) 



This stopping condition sets the value of V ma x. For the ra- 
dius R max , we adopt the geometric average of the two radii 
fulfilling equation Q. The origin of velocities is taken as 
the average velocity of all particles within 0.1-R m ax- Basic 
properties of our haloes are summarized in Table [T] 

The choice of -Rmax (and Vmax) has several advantages 
with respect to other alternatives, such as the virial radius or 
the radius r_2 at which the logarithmic slope of the density 
profile is equal to —2. It is well known that many proper- 
ties of dark matter haloes (e.g. the density profile) are not 
universal in terms of the virial radius only. Moreover, this 
quantity is not well defined for substructures, and tends to 
be unstable even for haloes undergoing weak interactions 
with close neighbours. While this does not pose a problem 
for the study of relaxed, isolated systems, it is not the opti- 
mal choice in the general case. The radius r_2, on the other 
hand, provides a much better reference point in the sense 
that many profiles display a universal structure when scaled 
in units of r_2, and it is much less sensitive to the environ- 
ment than the virial radius. However, it has the technical 
disadvantage that it is defined in terms of the local slope 
of the density profile, which is a second derivative of the 
mass distribution and therefore prone to significant numer- 
ical errors. This difficulty can be overcome by assuming an 
analytical model of the dark matter halo. In this case, R ma x 
and r_2 are proportional to one another, and the choice be- 
tween them becomes a matter of personal preference. We 
opted for R m ax because this quantity is model-independent; 
in fact, it does not rely on the assumption of a universal 
shape of the density profile. Moreover, the circular velocity 
of the halo is closer to an observable quantity than the halo 
mass, let alone the logarithmic slope of the density, and it 
reflects the depth of the potential well of the whole halo 
rather than a local property, which makes it a more robust 
estimator of the actual scale of the object. 



3 RESULTS 

In addition to the mass and circular velocity profiles, we 
obtained the dark matter density profile and its logarith- 
mic slope. The average radial velocity and the global angu- 
lar momentum are also computed within spherical shells, as 
well as the radial and tangential components of the velocity 
dispersion. This makes possible to investigate the radial de- 
pendence of other physical properties, such as the anisotropy 
parameter or the coarse-grained phase-space density. 

For all these quantities, we first consider the existence 
of a universal profile in terms of the two fundamental pa- 
rameters of the halo, which we take to be -R ma x and Vmax. 
When appropriate, we compare our numerical resul t s with 
the analytical m odels propos e d by iNavarro et al.l (|l997l . 
hereafter NFW), iRasia et all i|2004l . hereafter R04) and 
iDehnen fc McLaughlin! (|2005l . hereafter DM05) . In the NFW 
model, the maximum circular velocity occurs at -R max — 
2.163 r_2, while in R04, it takes place at -Rmax — 0.334 R v ir- 
For the DM05 model, we assume the quoted canonical val- 
ues, e = 3, /3o — —0.1 and Poo — 0.67. For this particular 
choice of the parameters, the maximum circular velocity oc- 
curs at a radius -Rmax — 1.852 ro, where ro denotes a char- 
acteristic radius of the halo. 

Possible systematic dependences on halo mass are in- 
vestigated by considering the results from our two sets of 
simulations (containing galaxy- and cluster-size objects at 
z = 0) separately. Temporal evolution is addressed by ana- 
lyzing the outputs at z = 0, 1, and 5. 



3.1 Mass distribution 

The circular velocity profiles of our 14 dark matter haloes, 
at the three different redshifts under study, are plotted to- 
gether on the left panel of Figure [T] Each individual profile 
is normalized at its maximum value, Vmax, an d the radius 
at which it is attained, R m ax- All haloes seem to be well 
described by a universal circular velocity profile; the scatter 
around the average, indicated by the error bars in Figure[T] is 
smaller than 30 per cent at all radii. This result, though not 
surprising, is remarkable in the sense that Vmax and R m ax 
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Figure 1. Circular velocity profiles, V c = \J GM/r, of our sample of dark matter haloes. Points with error bars show the average 
and one-cx scatter, including all objects, at all redshifts. On the left panel, individual profiles are plotted as green solid lines; dotted, 
dot-dashed, and dashed lines represent the NFW, R04 and DM05 models, respectively. On the middle and right panels, the global average 
(shown in green) is compared to the results obtained from our galaxy and cluster simulations. In both panels, solid, dashed and dotted 
lines correspond to the average profile at z = 0, 1, and 5, respectively. 




are not fits to the dark matter distribution, but physical 
properties of the halo. Indeed, the radius i? max (as well as 
r-2) corresponds to the transition between two completely 
different dynamical regime^], and, particularly at high red- 
shift, we claim it is a much better indicator of the region 
within which the halo can be considered in equilibrium. A 
first hint into that direction can already be seen on the mid- 
dle and right panels, which show the evolution of galaxies 
and (proto)clusters, respectively. Whereas cosmic evolution 
plays an important role outside -R max , the average behavior 
of the circular velocity profile is much more homogeneous 
within this radius. 

A similar conclusion can be reached in terms of the den- 
sity profile. The density profiles of our dark matter haloes 
are plotted in Figure [2] normalized to the characteristic den- 
sity 



1 For r < -Rmax, the mass is dominated by particles with apoc- 
entre a > r (i.e. coming from out er regions). Outside fl max , 
most particles hav e a < r (see e.g. iFillmore fc GoIdreicb|[l984l ; 
lBertschinger1ll985l) . 



pmax — 




and the radius -Rmax. Again, neither the density profile nor 
its logarithmic slope 7 (plotted on the top sub-panels) devi- 
ate significantly from the average behaviour for r < J? max . 
In the outer regions, one does not only observe the effects of 
cosmic evolution, but also a significant variation from object 
to object (much larger than any dependence on halo mass) 
which reflects the different environments in which our dark 
matter haloes are embedded. Within i? max , we do not ob- 
serve any systematic trend with mass or redshift, but one 
should bear in mind that small departures from universal- 
ity (in particular, smaller than the scatter) could only be 
detected with a much larger sample of simulated objects. 

Concerning the precise shape of the density or circular 
velocity profiles, our results are consistent up to 7? max with 
all of the three analytical models considered. The average 
density profile of our haloes seems to be somewhat steeper 
than the NFW formula, and it is better reproduced by the 
functions proposed by R04 or DM05. This, however, is prob- 
ably related to the fact that we have not made any selection 
based on the dynamical state of our haloes, while the NFW 
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model is intended to provide a good description of relaxed 
objects (see e.g. the discussion in R04). 



3.2 Bulk motion 

The average velocity field of dark matter haloes can be sep- 
arated into a radial and a tangential component; radial mo- 
tion tells us about the infall of matter into the halo, while 
tangential motion, measured by the spin parameter, contains 
information about the recent accretion history and the sur- 
rounding tidal field. On the other hand, the offsets between 
the coordinates (position and velocity) of the centre of mass 
on different scales reflect the displacement of the halo with 
respect to the nearby large-scale structures. 

Let us consider first the motion of the dark matter halo 
as a whole. If haloes were smooth, spherically-symmetric sys- 
tems in perfect equilibrium, the position of centre of mass 
should be independent of the scale considered. This state- 
ment would also be true in triaxial, rotating systems, and it 
holds even in the case of significant infall, as long as the mass 
contributed by substructures is small compared to the total 
mass of the material being accreted. Therefore, one would 
expect that a possible offset between the central and outer 
regions of the halo should be small in relaxed systems. Along 
these lines, the distance between the most bound particle 
(close to our adopted centre) and the centre of mass evalu- 



ated at the scale of the virial radius has been used to detect 
unrelaxe d haloes and quantify departures fro m equilibrium 
(see e.g. iMaccio et ai1l2007l ; iNeto et alj|2007l ). 

Figure[3]shows the offset of the position and the velocity 
of the centre of mass (evaluated within a sphere of radius 
r) with respect to our adopted origin of coordinates, which 
traces the highest density peak. The behaviour of the centre 
of mass velocity (plotted on the bottom panels) is easy to 
understand: the offset is minimum at r = 0.1J? m ax because 
this is exactly the prescription we followed to set up the 
origin (see Section [2}. At smaller radii, we detect the motion 
of the very central regions, which is typically of the order of 
ten per cent of the maximum circular velocity. At large radii, 
the velocity of the centre of mass tends asymptotically to the 
average velocity of the halo with respect to the background 
universe. The same pattern is observed for our galaxy and 
cluster simulations, at the three redshifts considered. 

The behaviour of the centre of mass position is at first 
sight more striking: rather than converging to an asymp- 
totic value, the offset increases roughly linearly with radius, 
up to scales that can be as large as the box size for the 
largest systems. Again, the same pattern is common to all 
systems, at all epochs, although the change in normaliza- 
tion with redshift seems to indicate that -R max is probably 
not the relevant scaling radius. Obviously, the offset in po- 
sition should stay constant beyond the scale of homogene- 
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Figure 6. Radial velocity dispersion. Line styles are the same as in Figure IT1 



ity, but this is larger than the 80 Mpc of the simulation 
box. What Figure [3] shows is that, on the scales probed by 
our numerical experiments, the universe is not entirely ho- 
mogenous, and the tid al field of the large-scale struc ture is 
not negligible (see e.g. iMartinez-Vaquero et al.l [2007) . One 
should bear in mind, however, that the finite box size and 
periodic boundary conditions severely limit our ability to 
sample the Fourier modes of the primordial fluctuations on 
large scales. For L ~ 80 h~ Mpc, the amplitude of the CDM 
power spectrum is not negligible, and therefore it is possible 
that our results on such scales are not representative of the 
real universe. 

Radial velocity profiles are plotted in Figure [4] Al- 
though there is, in an approximate sense, a general pattern 
common to all objects, at all epochs (an inner region in equi- 
librium, an infall region, and a region still expanding), the 
details depend a lot on the scale of the halo and the redshift 
of observation. In units of i? max , the turn-around radius of 
all systems stays approximatel y constant a t J? ta ~ 5i? max , 
but we confirm the claims by Prada et al.l (|2006h that, at 
2 = 0, the infall region is effectively absent in galaxy-mass 
haloes. At earlier times (z = 1), the progenitors of galaxies 
display an infall pattern similar to that of clusters at z — 0, 
while the progenitors of present-day clusters feature much 
higher infall velocities. In both cases, the equilibrium region 
where v T w is closer to i? max than to the virial radius. The 
importance of infall becomes even more dramatic at z = 5. 



Maximum infall velocities are found at r ~ 2i? max (about or 
within the virial radius of the halo), with magnitudes com- 
parable to the maximum circular velocity of the object. The 
transition to equilibrium takes place at r ~ i? max for proto- 
galaxies and r ~ 0.3i? max for proto-clusters, relatively close 
to our resolution limit. 

We measure the total angular momentu m of the halo 
in ter ms of the spin parameter, defined as in iBullock et al.l 



A(r) = 



J 



V2GM 3 



1 «tan 

y/2 Vc 



(6) 



where wtan = J /Mr is the average tangential velocity within 
radius r. The distribution of A in our 14 haloes is con- 
sistent with the average value Ao ~ 0.035 and the spread 
A ln(A) ~ 0.6 that are usually reported at the virial radius. 
Indeed, we find (see Figure [5]) that such values are roughly 
representative of the bulk rotation of the halo at all radii, al- 
beit there may be large deviations in the individual profiles. 
Nonetheless, we also observe a systematic trend in most ob- 
jects, as well as on the average profiles: the spin parameter 
increases with radius until it reaches a maximum value, and 
then it decreases again towards the outer parts. The inten- 
sity and location of this maximum are strongly dependent 
on epoch, and perhaps on the mass of the object, the maxi- 
mum being stronger and closer to the centre as we consider 
more massive systems at higher redshift. 
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Figure 7. Anisotropy parameter (top) and and its relation with the local logarithmic slope of the density profile (bottom). Line styles 
are the same as in Figure [T] Open squares on the bottom left panel represent individual measurements of 7 and at a given radius 
(only data whithin iJ ma x are considered), and vertical lines indicate the range used to fit the f} — 7 relation. 



3.3 Random motion 

In general, the average bulk motion of the halo is much 
slower than the velocity of individual dark matter particles. 
The latter, which determines the orbital structure and there- 
fore the dynamics of the system, is dominated by random 
motions. 

Radial velocity dispersion profiles are plotted in Fig- 
ure [6] Again, we observe a remarkable degree of universality 
in terms of i? max and Knax, with a scatter of about 20 per 
cent with respect to the average behaviour within -R max . Nei- 
ther the scale of the object nor the cosmic epoch seem to 
have a significant impact on the results, except for a slight 
steepening outside -Rmax at high redshift. The radial velocity 
dispersion is well described by the R04 and DM05 models 
in the inner regions, but both of them seem to be system- 
atically higher than our measurements in the outskirts of 
haloes (near the virial radius). We note, however, that the 
discrepancy is only of the order of one sigma, and therefore 
more statistics are required to verify whether it is real or 
not. 

The anisotropy parameter, defined as 

/3 = l-gf (7) 

relates the radial and tangential components of the velocity 
dispersion. For an isotropic system, (3 — 0, while /3 = 1 and 
P = — 00 correspond to purely radial and circular orbits, 
respectively. Our results, shown on the top panels of Fig- 



ure [7] are well reproduced by both R04 and DM05 models 
within -Rmax- Most of the large scatter observed arises from 
noise in the numerical evaluation of /3, but part of it (espe- 
cially in the outer regions) is also due to the environment 
of individual haloes. For all object types except galaxies at 
2 = (the only ones that do not display an infall region), 
the anisotropy parameter attains a maximum around i? max , 
and the decline towards larger radii seems to be sharper at 

high redshift than at the present t ime. 

It has been su ggested (|Hansen fc Moorel [2006; 
lHansen fc Stadell [2006) that a linear relation exists 
between the anisotropy parameter and the local value 
of the logarithmic slope of the density profile, 7. As can 
be seen on the bottom panels of Figure [7] our data are 
compatible with such a relation. A least-squares fit of the 
form 

P = A + B 7 (8) 

over the range -2.5 < 7 < -1 yields A = 0.016 ± 0.016 
and B — —0.085 ± 0.009, with a correlation coefficient r = 
—0.97. In this region, such a fit does not deviate significantly 
from the predictions of the R04 a nd DM05 models nor th e 
results of lHansen fc Moorel i|2006l ); lHansen fc Stadel (|2006h . 
although the slope and intercept we obtain are different from 
the values quoted by these authors, A = —0.16 and B — 
—0.2. Obviously, the extrapolation outside the fitted range 
(which we would strongly discourage) leads to substantially 
different results in both directions. 



s 
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Finally, we investigate the coarse-grained phase-space 
density profile of our haloes in Figure [8] comparing three 
different definitions of this quantity. On the top panels, we 
plot 

<7tot (r) = -=- 



(9) 

(3 — 2/9)ctj denotes the total veloc- 
ity dispersion at radius r. On the middle panels, the radial 
velocity dispersion is used instead, 



where ov 



„2 | J2 



qr{r) 



(3 - 2(3) 3/2 g tot 



(10) 



and on the bottom panel, the coarse-grained phase-space 
density is defined in terms of cumulative quantities, 



«~(r) = 4- 



(11) 



where (p) = 3M/inr 3 is the average density within radius 
r and of os = J (v 2 — {v) 2 )^nr 2 p dr is the average line-of- 
sight velocity dispersion. 

Our results indicate that the coarse-grained phase-space 
density profile, defined in any of the aforementioned ways, 
can be approximately described by a power law, in agree- 
ment with previous work. Considering data in the range 
0.1 < r/Rmax < 1, we find 

qi(r) 



Qn 



= K ir ~ 



(12) 



with best-fitting exponents Qtot = 1-77 ± 0.04, ct T — 1.87 ± 
0.03, and Qi os = 1.82±0.01 for the quantities gtot, Qi an d q\ os , 
respectively. The corresponding normalizations are K tot = 

10 -0.50±0.0 3) ^ = yjO.OTiO.CB ^ ios = 10 O.51±O.Ol 5 and 

the quantity Q ma x is defined as 

Pmax 3 



Q 



(13) 



in all cases. A qualitatively similar behaviour is obtained 
for p/aLn, with a tan = 1.70±0.05 and iC tan = io-° 15±0 03 . 

Regardless of the exact definition of q, both sets of sim- 
ulations (galaxy- and cluster-size objects) follow the same 
power law, and cosmic evolution can only be appreciated 
at very large radii (usually well beyond the virial radius jfj. 
Therefore, these power laws can be used to provide a simple, 
phenomenological description of the structure of dark mat- 
ter haloes. However, as in the case of the anisotropy profile 
(or any of the other relations investigated here) , we would 
strongly advise against assuming that any given analytic fit 
is an exact description, or that it can be safely extrapolated 
outside its range of validity. 

For instance, the first two definitions of the coarse- 
grained phase space density could only follow exactly a 
power law if the anisotropy was constant at all radii, which 
is definitely not the case, or it varied according to 



0(r) = 



(r/RiY 



(14) 



where the exponent e = a T — a to t — 0.1 would depend on 
the precise values of the logarithmic slopes of the q T and gtot 
profiles, and Ri would be a characteristic radius at which 
f3 = 1 and q r = q to t- Expression (|14[) does indeed provide a 
good fit to our simulated data, but one should not conclude 
from there that the innermost orbits are perfectly circular, 
i.e. (5(0) = — oo, because there is no a priori reason to think 
that q r (r) and gtot(^) ought to be exact power laws. 



2 The apparent deviations at r < 0.01i? m ax should be taken 
with a grain of salt, since they might well be related to our choice 



of the origin of velocities. 
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4 CONCLUSIONS 

In this paper, we have studied the internal structure of 
dark matter haloes using high-resolution cosmological N- 
body simulations of individual objects. Our sample consists 
of six dark matter haloes on galactic scales and eight systems 
that correspond to galaxy groups and clusters at the present 
day. These objects have been randomly selected, without 
imposing any condition on their dynamical state. We have 
considered the radially-averaged profiles of cumulative mass 
and circular velocity, dark matter density and its logarith- 
mic slope, position and velocity offsets of the centre of mass 
on different scales, average radial velocity and spin parame- 
ter, radial velocity dispersion and anisotropy parameter, and 
coarse-grained phase-space density. 

Our results show that the all these profiles display a 
similar structure when expressed in terms of the maximum 
circular velocity of the halo, Knax, and the radius at which 
it is attained, 7? max . Within -R max , this structure is inde- 
pendent on object scale and cosmic epoch, and it is well 
fit by several analytical models that have previously been 
proposed to describe the numerical data at z — in terms 
of other scaling parameters. This indicates that the equi- 
librium properties of dark matter haloes are independent of 
the individual merger histories, and therefore they must be 
set by some physical process or processes (such as violent re- 
laxation, phase mixing, radial orbit instability, etc.) whose 
outcome is similar for all objects. In particular, the present 
study rules out theoretical models that predict mass or red- 
shift dependence of the dynamical structure of dark matter 
haloes. 

The radius 7? max gives a more conservative indication 
of the region where a given object can be considered in equi- 
librium than other measurements, such as the virial radius. 
This is particularly important at high redshift, where sig- 
nificant radial infall (up to fifty per cent of Knax) may take 
place within the latter. Outside f? max , not only the radial 
velocity profile, but also the mass distribution, the spin, or 
the anisotropy parameter display strong systematic depen- 
dences on mass and/or redshift. It is in this outer region 
where the system retains some memory of the primordial 
initial conditions and displays more obvious signatures of 
its current environment. 

Finally, we provide fits for two relations between 
position- and velocity-related quantities. We corroborate the 
roughly linear relation between the anisotropy parameter 
and the local logarithmic slope of the density profile, al- 
though the best-fitting zero-point and proportionality con- 
stant we obtain differ from the values originally proposed. 
We also compute the slope and normalization of the coarse- 
grained phase-space density profile, which is found to be 
fairly well described by a pure power law, according to 
three different definitions. These two quantities, as well 
as all the others considered in the present study, involve 
integrals of the underlying distribution function (the fine- 
grained, six-dimensional phase-space density) over the ve- 
locity. The distribut ion function itse l f is very difficult to 



measure accurately (lArad et al.1 1 20041 ; lAscasibar fc Binnevl 
l2005l ; ISharma fc Steinmet j|2006f) . 



but it should certainly be 
the aim of future work in order to fully characterize the 
structure of dark matter haloes. Promising steps in th at di- 
rection have been taken by e.g. lHansen et al.f (|2006l ). who 



consider the velocity distrib ution function withi n different 
bins in potential energy, and lWoitak et all (|2008T l. who pro- 
pose a simple parametric form in terms of the energy, angu- 
lar momentum, and anisotropy of particle orbits. 

We would like to stress, though, that caution must be 
exercised when interpreting the results of phenomenological 
studies (like the present one). In particular, analytical fits 
are an extremely useful tool that provides an approximate 
description of dark matter haloes, but they constitute by no 
means a physical law. Extrapolating outside the appropriate 
validity range of or using a fit as an exact description can 
lead to misleading conclusions, not only quantitatively, but 
also on a qualitative level. More generally, we do not share 
the widespread opinion that the formation and evolution of 
dark matter haloes is nowadays well understood. Even if we 
could provide a good phenomenological prescription of their 
dynamical structure in terms of empirical relations, there is 
an absolute lack of understanding of their physical origin. 
Open questions are, for example, how many parameters are 
actually needed to describe a dark matter halo? Which re- 
lations (if any) are indeed fundamental? What physical pro- 
cesses drive the evolution of dark matter haloes, and why do 
they lead to a 'universal' equilibrium state? We hope that 
the results presented here help in providing a benchmark to 
test any physically-motivated model that is able to address 
these questions. 
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